! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! The ODE Function of Chemical Model File
! 
! Generated by KPP-2.2.3 symbolic chemistry Kinetics PreProcessor
!       (http://www.cs.vt.edu/~asandu/Software/KPP)
! KPP is distributed under GPL, the general public licence
!       (http://www.gnu.org/copyleft/gpl.html)
! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa
! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech
!     With important contributions from:
!        M. Damian, Villanova University, USA
!        R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany
! 
! File                 : aromatics_kpp_Function.f90
! Time                 : Mon Nov 23 18:26:23 2020
! Working directory    : /n/home08/kbates/Aromatics/CRI_chamber
! Equation file        : aromatics_kpp.kpp
! Output root filename : aromatics_kpp
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



MODULE aromatics_kpp_Function

  USE aromatics_kpp_Parameters
  IMPLICIT NONE

! A - Rate for each equation
  REAL(kind=dp) :: A(NREACT)

CONTAINS


! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
! 
! Fun - time derivatives of variables - Agregate form
!   Arguments :
!      V         - Concentrations of variable species (local)
!      F         - Concentrations of fixed species (local)
!      RCT       - Rate constants (local)
!      Vdot      - Time derivative of variable species concentrations
! 
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Fun ( V, F, RCT, Vdot )

! V - Concentrations of variable species (local)
  REAL(kind=dp) :: V(NVAR)
! F - Concentrations of fixed species (local)
  REAL(kind=dp) :: F(NFIX)
! RCT - Rate constants (local)
  REAL(kind=dp) :: RCT(NREACT)
! Vdot - Time derivative of variable species concentrations
  REAL(kind=dp) :: Vdot(NVAR)


! Computation of equation rates
  A(1) = RCT(1)*V(78)
  A(2) = RCT(2)*V(78)*V(193)
  A(3) = RCT(3)*V(78)*V(189)
  A(4) = RCT(4)*V(78)*V(192)
  A(5) = RCT(5)*V(78)*V(192)
  A(6) = RCT(6)*V(19)
  A(7) = RCT(7)*V(189)*V(193)
  A(8) = RCT(8)*V(192)*V(193)
  A(9) = RCT(9)*V(189)*V(189)
  A(10) = RCT(10)*V(189)*V(195)
  A(11) = RCT(11)*V(192)*V(195)
  A(12) = RCT(12)*V(192)*V(195)
  A(13) = 5.35e+07*V(19)
  A(14) = RCT(14)*V(193)*V(197)
  A(15) = RCT(15)*V(56)*V(197)
  A(16) = RCT(16)*V(82)*V(197)
  A(17) = RCT(17)*V(193)*V(196)
  A(18) = RCT(18)*V(196)*V(197)
  A(19) = RCT(19)*V(196)*V(196)
  A(20) = RCT(20)*V(189)*V(197)
  A(21) = RCT(21)*V(192)*V(197)
  A(22) = 2e-11*V(195)*V(197)
  A(23) = RCT(23)*V(189)*V(196)
  A(24) = RCT(24)*V(192)*V(196)
  A(25) = RCT(25)*V(63)*V(197)
  A(26) = 4e-12*V(195)*V(196)
  A(27) = RCT(27)*V(28)*V(197)
  A(28) = RCT(28)*V(145)*V(197)
  A(29) = RCT(29)*V(193)
  A(30) = RCT(30)*V(193)
  A(31) = RCT(31)*V(82)
  A(32) = RCT(32)*V(192)
  A(33) = RCT(33)*V(195)
  A(34) = RCT(34)*V(195)
  A(35) = RCT(35)*V(28)
  A(36) = RCT(36)*V(145)
  A(37) = RCT(37)*V(195)*V(195)
  A(38) = RCT(38)*V(22)
  A(39) = RCT(39)*V(63)
  A(40) = RCT(40)*V(22)
  A(41) = RCT(41)*V(63)
  A(42) = RCT(42)*V(30)*V(197)
  A(43) = RCT(43)*V(189)*V(190)
  A(44) = RCT(44)*V(190)*V(196)
  A(45) = RCT(45)*V(190)*V(190)
  A(46) = RCT(46)*V(190)*V(190)
  A(47) = RCT(47)*V(190)*V(192)
  A(48) = RCT(48)*V(50)
  A(49) = RCT(49)*V(52)*V(197)
  A(50) = RCT(50)*V(52)*V(197)
  A(51) = RCT(51)*V(174)*V(197)
  A(52) = RCT(52)*V(174)*V(195)
  A(53) = RCT(53)*V(152)*V(197)
  A(54) = RCT(54)*V(194)*V(197)
  A(55) = RCT(55)*V(52)
  A(56) = RCT(56)*V(174)
  A(57) = RCT(57)*V(174)
  A(58) = RCT(58)*V(50)
  A(59) = RCT(59)*V(50)
  A(60) = RCT(60)*V(90)*V(197)
  A(61) = RCT(61)*V(90)*V(197)
  A(62) = RCT(62)*V(96)*V(197)
  A(63) = RCT(63)*V(177)*V(197)
  A(64) = RCT(64)*V(177)*V(195)
  A(65) = RCT(65)*V(191)*V(192)
  A(66) = RCT(66)*V(76)
  A(67) = RCT(67)*V(189)*V(191)
  A(68) = RCT(68)*V(41)*V(197)
  A(69) = RCT(69)*V(170)*V(189)
  A(70) = RCT(70)*V(20)*V(197)
  A(71) = RCT(71)*V(20)*V(197)
  A(72) = RCT(72)*V(144)*V(189)
  A(73) = RCT(73)*V(159)*V(189)
  A(74) = RCT(74)*V(40)*V(197)
  A(75) = RCT(75)*V(183)*V(189)
  A(76) = RCT(76)*V(183)*V(189)
  A(77) = RCT(77)*V(172)*V(189)
  A(78) = RCT(78)*V(166)*V(189)
  A(79) = RCT(79)*V(167)*V(189)
  A(80) = RCT(80)*V(184)*V(189)
  A(81) = RCT(81)*V(184)
  A(82) = RCT(82)*V(169)*V(189)
  A(83) = RCT(83)*V(179)*V(189)
  A(84) = RCT(84)*V(179)
  A(85) = RCT(85)*V(161)*V(189)
  A(86) = RCT(86)*V(154)*V(189)
  A(87) = RCT(87)*V(181)*V(189)
  A(88) = RCT(88)*V(181)*V(195)
  A(89) = RCT(89)*V(176)*V(195)
  A(90) = RCT(90)*V(163)*V(195)
  A(91) = RCT(91)*V(163)*V(189)
  A(92) = RCT(92)*V(163)*V(192)
  A(93) = RCT(93)*V(163)*V(196)
  A(94) = RCT(94)*V(153)*V(195)
  A(95) = RCT(95)*V(153)*V(189)
  A(96) = RCT(96)*V(153)*V(196)
  A(97) = RCT(97)*V(153)*V(190)
  A(98) = RCT(98)*V(176)*V(193)
  A(99) = RCT(99)*V(176)*V(197)
  A(100) = RCT(100)*V(148)*V(189)
  A(101) = RCT(101)*V(148)*V(196)
  A(102) = RCT(102)*V(148)*V(190)
  A(103) = RCT(103)*V(181)*V(181)
  A(104) = RCT(104)*V(156)*V(189)
  A(105) = RCT(105)*V(40)*V(195)
  A(106) = RCT(106)*V(185)*V(197)
  A(107) = RCT(107)*V(137)*V(197)
  A(108) = RCT(108)*V(187)*V(197)
  A(109) = RCT(109)*V(171)*V(192)
  A(110) = RCT(110)*V(24)
  A(111) = RCT(111)*V(51)*V(192)
  A(112) = RCT(112)*V(175)*V(192)
  A(113) = RCT(113)*V(129)
  A(114) = RCT(114)*V(113)
  A(115) = RCT(115)*V(171)*V(189)
  A(116) = RCT(116)*V(175)*V(189)
  A(117) = RCT(117)*V(187)*V(195)
  A(118) = RCT(118)*V(115)*V(197)
  A(119) = RCT(119)*V(144)*V(190)
  A(120) = RCT(120)*V(159)*V(190)
  A(121) = RCT(121)*V(183)*V(196)
  A(122) = RCT(122)*V(172)*V(196)
  A(123) = RCT(123)*V(166)*V(196)
  A(124) = RCT(124)*V(167)*V(196)
  A(125) = RCT(125)*V(184)*V(196)
  A(126) = RCT(126)*V(169)*V(196)
  A(127) = RCT(127)*V(179)*V(196)
  A(128) = RCT(128)*V(161)*V(196)
  A(129) = RCT(129)*V(154)*V(196)
  A(130) = RCT(130)*V(181)*V(196)
  A(131) = RCT(131)*V(156)*V(196)
  A(132) = RCT(132)*V(188)*V(197)
  A(133) = RCT(133)*V(170)*V(190)
  A(134) = RCT(134)*V(188)*V(195)
  A(135) = RCT(135)*V(183)*V(190)
  A(136) = RCT(136)*V(172)*V(190)
  A(137) = RCT(137)*V(166)*V(190)
  A(138) = RCT(138)*V(167)*V(190)
  A(139) = RCT(139)*V(184)*V(190)
  A(140) = RCT(140)*V(184)*V(184)
  A(141) = RCT(141)*V(173)*V(190)
  A(142) = RCT(142)*V(169)*V(190)
  A(143) = RCT(143)*V(179)*V(190)
  A(144) = RCT(144)*V(161)*V(190)
  A(145) = RCT(145)*V(154)*V(190)
  A(146) = RCT(146)*V(181)*V(190)
  A(147) = RCT(147)*V(156)*V(190)
  A(148) = RCT(148)*V(66)*V(197)
  A(149) = RCT(149)*V(120)*V(197)
  A(150) = RCT(150)*V(170)*V(170)
  A(151) = RCT(151)*V(170)*V(170)
  A(152) = RCT(152)*V(170)*V(196)
  A(153) = RCT(153)*V(144)*V(196)
  A(154) = RCT(154)*V(159)*V(196)
  A(155) = RCT(155)*V(191)*V(196)
  A(156) = RCT(156)*V(171)*V(196)
  A(157) = RCT(157)*V(175)*V(196)
  A(158) = RCT(158)*V(143)*V(197)
  A(159) = RCT(159)*V(143)*V(193)
  A(160) = RCT(160)*V(113)*V(197)
  A(161) = RCT(161)*V(129)*V(197)
  A(162) = RCT(162)*V(113)*V(193)
  A(163) = RCT(163)*V(1)*V(197)
  A(164) = RCT(164)*V(165)*V(197)
  A(165) = RCT(165)*V(165)*V(197)
  A(166) = RCT(166)*V(143)*V(195)
  A(167) = RCT(167)*V(155)*V(197)
  A(168) = RCT(168)*V(178)*V(197)
  A(169) = RCT(169)*V(155)*V(195)
  A(170) = RCT(170)*V(178)*V(195)
  A(171) = RCT(171)*V(142)*V(197)
  A(172) = RCT(172)*V(197)*V(198)
  A(173) = RCT(173)*V(186)*V(197)
  A(174) = RCT(174)*V(182)*V(197)
  A(175) = RCT(175)*V(182)*V(197)
  A(176) = RCT(176)*V(144)*V(191)
  A(177) = RCT(177)*V(159)*V(191)
  A(178) = RCT(178)*V(144)*V(191)
  A(179) = RCT(179)*V(159)*V(191)
  A(180) = RCT(180)*V(142)*V(193)
  A(181) = RCT(181)*V(193)*V(198)
  A(182) = RCT(182)*V(186)*V(193)
  A(183) = RCT(183)*V(142)*V(195)
  A(184) = RCT(184)*V(186)*V(195)
  A(185) = RCT(185)*V(186)*V(195)
  A(186) = RCT(186)*V(171)*V(190)
  A(187) = RCT(187)*V(175)*V(190)
  A(188) = RCT(188)*V(171)*V(190)
  A(189) = RCT(189)*V(175)*V(190)
  A(190) = RCT(190)*V(112)*V(197)
  A(191) = RCT(191)*V(112)*V(197)
  A(192) = RCT(192)*V(81)*V(197)
  A(193) = RCT(193)*V(53)*V(197)
  A(194) = RCT(194)*V(54)*V(197)
  A(195) = RCT(195)*V(55)*V(197)
  A(196) = RCT(196)*V(84)*V(197)
  A(197) = RCT(197)*V(77)*V(197)
  A(198) = RCT(198)*V(80)*V(197)
  A(199) = RCT(199)*V(102)*V(197)
  A(200) = RCT(200)*V(92)*V(197)
  A(201) = RCT(201)*V(103)*V(197)
  A(202) = RCT(202)*V(102)*V(197)
  A(203) = RCT(203)*V(92)*V(197)
  A(204) = RCT(204)*V(103)*V(197)
  A(205) = RCT(205)*V(2)*V(197)
  A(206) = RCT(206)*V(31)*V(197)
  A(207) = RCT(207)*V(32)*V(197)
  A(208) = RCT(208)*V(29)*V(197)
  A(209) = RCT(209)*V(151)*V(196)
  A(210) = RCT(210)*V(151)*V(189)
  A(211) = RCT(211)*V(95)*V(197)
  A(212) = RCT(212)*V(93)*V(197)
  A(213) = RCT(213)*V(97)*V(197)
  A(214) = RCT(214)*V(97)*V(197)
  A(215) = RCT(215)*V(108)*V(197)
  A(216) = RCT(216)*V(108)*V(197)
  A(217) = RCT(217)*V(164)*V(191)
  A(218) = RCT(218)*V(164)*V(191)
  A(219) = RCT(219)*V(164)*V(190)
  A(220) = RCT(220)*V(164)*V(164)
  A(221) = RCT(221)*V(164)*V(196)
  A(222) = RCT(222)*V(164)*V(189)
  A(223) = RCT(223)*V(164)*V(189)
  A(224) = RCT(224)*V(42)*V(197)
  A(225) = RCT(225)*V(41)*V(195)
  A(226) = RCT(226)*V(180)*V(193)
  A(227) = RCT(227)*V(191)*V(191)
  A(228) = RCT(228)*V(190)*V(191)
  A(229) = RCT(229)*V(190)*V(191)
  A(230) = RCT(230)*V(183)*V(191)
  A(231) = RCT(231)*V(166)*V(191)
  A(232) = RCT(232)*V(167)*V(191)
  A(233) = RCT(233)*V(184)*V(191)
  A(234) = RCT(234)*V(173)*V(191)
  A(235) = RCT(235)*V(169)*V(191)
  A(236) = RCT(236)*V(179)*V(191)
  A(237) = RCT(237)*V(154)*V(191)
  A(238) = RCT(238)*V(172)*V(191)
  A(239) = RCT(239)*V(161)*V(191)
  A(240) = RCT(240)*V(181)*V(191)
  A(241) = RCT(241)*V(156)*V(191)
  A(242) = RCT(242)*V(183)*V(191)
  A(243) = RCT(243)*V(166)*V(191)
  A(244) = RCT(244)*V(167)*V(191)
  A(245) = RCT(245)*V(184)*V(191)
  A(246) = RCT(246)*V(173)*V(191)
  A(247) = RCT(247)*V(169)*V(191)
  A(248) = RCT(248)*V(179)*V(191)
  A(249) = RCT(249)*V(172)*V(191)
  A(250) = RCT(250)*V(161)*V(191)
  A(251) = RCT(251)*V(181)*V(191)
  A(252) = RCT(252)*V(156)*V(191)
  A(253) = RCT(253)*V(154)*V(191)
  A(254) = RCT(254)*V(170)*V(191)
  A(255) = RCT(255)*V(170)*V(191)
  A(256) = RCT(256)*V(171)*V(191)
  A(257) = RCT(257)*V(175)*V(191)
  A(258) = RCT(258)*V(134)*V(193)
  A(259) = RCT(259)*V(168)*V(193)
  A(260) = RCT(260)*V(180)*V(197)
  A(261) = RCT(261)*V(173)*V(189)
  A(262) = RCT(262)*V(173)*V(189)
  A(263) = RCT(263)*V(173)*V(196)
  A(264) = RCT(264)*V(168)*V(197)
  A(265) = RCT(265)*V(134)*V(197)
  A(266) = RCT(266)*V(157)*V(189)
  A(267) = RCT(267)*V(157)*V(196)
  A(268) = RCT(268)*V(162)*V(189)
  A(269) = RCT(269)*V(162)*V(196)
  A(270) = RCT(270)*V(104)*V(197)
  A(271) = RCT(271)*V(147)*V(197)
  A(272) = RCT(272)*V(158)*V(197)
  A(273) = RCT(273)*V(126)*V(189)
  A(274) = RCT(274)*V(126)*V(196)
  A(275) = RCT(275)*V(126)*V(192)
  A(276) = RCT(276)*V(23)
  A(277) = RCT(277)*V(101)*V(197)
  A(278) = RCT(278)*V(123)*V(189)
  A(279) = RCT(279)*V(123)*V(189)
  A(280) = RCT(280)*V(123)*V(196)
  A(281) = RCT(281)*V(100)*V(197)
  A(282) = RCT(282)*V(122)*V(189)
  A(283) = RCT(283)*V(122)*V(189)
  A(284) = RCT(284)*V(122)*V(196)
  A(285) = RCT(285)*V(100)*V(193)
  A(286) = RCT(286)*V(109)*V(197)
  A(287) = RCT(287)*V(160)*V(197)
  A(288) = RCT(288)*V(141)*V(194)
  A(289) = RCT(289)*V(141)*V(189)
  A(290) = RCT(290)*V(141)*V(192)
  A(291) = RCT(291)*V(141)*F(1)
  A(292) = RCT(292)*V(141)*F(1)
  A(293) = RCT(293)*V(133)*V(194)
  A(294) = RCT(294)*V(133)*V(189)
  A(295) = RCT(295)*V(133)*V(192)
  A(296) = RCT(296)*V(133)*F(1)
  A(297) = RCT(297)*V(133)*F(1)
  A(298) = RCT(298)*V(110)*V(194)
  A(299) = RCT(299)*V(110)*V(189)
  A(300) = RCT(300)*V(110)*V(192)
  A(301) = RCT(301)*V(110)*F(1)
  A(302) = RCT(302)*V(121)*V(194)
  A(303) = RCT(303)*V(121)*V(189)
  A(304) = RCT(304)*V(121)*V(192)
  A(305) = RCT(305)*V(121)*F(1)
  A(306) = RCT(306)*V(127)*V(194)
  A(307) = RCT(307)*V(127)*V(189)
  A(308) = RCT(308)*V(127)*V(192)
  A(309) = RCT(309)*V(127)*F(1)
  A(310) = RCT(310)*V(127)*F(1)
  A(311) = RCT(311)*V(119)*V(194)
  A(312) = RCT(312)*V(119)*V(189)
  A(313) = RCT(313)*V(119)*V(192)
  A(314) = RCT(314)*V(119)*F(1)
  A(315) = RCT(315)*V(119)*F(1)
  A(316) = RCT(316)*V(114)*V(194)
  A(317) = RCT(317)*V(114)*V(189)
  A(318) = RCT(318)*V(114)*V(192)
  A(319) = RCT(319)*V(114)*F(1)
  A(320) = RCT(320)*V(114)*F(1)
  A(321) = RCT(321)*V(98)*V(197)
  A(322) = RCT(322)*V(87)*V(197)
  A(323) = RCT(323)*V(91)*V(189)
  A(324) = RCT(324)*V(94)*V(189)
  A(325) = RCT(325)*V(91)*V(196)
  A(326) = RCT(326)*V(94)*V(196)
  A(327) = RCT(327)*V(91)
  A(328) = RCT(328)*V(91)
  A(329) = RCT(329)*V(177)
  A(330) = RCT(330)*V(76)
  A(331) = RCT(331)*V(187)
  A(332) = RCT(332)*V(115)
  A(333) = RCT(333)*V(115)
  A(334) = RCT(334)*V(188)
  A(335) = RCT(335)*V(165)
  A(336) = RCT(336)*V(155)
  A(337) = RCT(337)*V(155)
  A(338) = RCT(338)*V(155)
  A(339) = RCT(339)*V(178)
  A(340) = RCT(340)*V(198)
  A(341) = RCT(341)*V(198)
  A(342) = RCT(342)*V(198)
  A(343) = RCT(343)*V(186)
  A(344) = RCT(344)*V(182)
  A(345) = RCT(345)*V(112)
  A(346) = RCT(346)*V(81)
  A(347) = RCT(347)*V(53)
  A(348) = RCT(348)*V(54)
  A(349) = RCT(349)*V(55)
  A(350) = RCT(350)*V(84)
  A(351) = RCT(351)*V(80)
  A(352) = RCT(352)*V(77)
  A(353) = RCT(353)*V(102)
  A(354) = RCT(354)*V(92)
  A(355) = RCT(355)*V(103)
  A(356) = RCT(356)*V(95)
  A(357) = RCT(357)*V(104)
  A(358) = RCT(358)*V(93)
  A(359) = RCT(359)*V(97)
  A(360) = RCT(360)*V(108)
  A(361) = RCT(361)*V(185)
  A(362) = RCT(362)*V(42)
  A(363) = RCT(363)*V(158)
  A(364) = RCT(364)*V(147)
  A(365) = RCT(365)*V(134)
  A(366) = RCT(366)*V(168)
  A(367) = RCT(367)*V(160)
  A(368) = RCT(368)*V(90)
  A(369) = RCT(369)*V(129)
  A(370) = RCT(370)*V(109)
  A(371) = RCT(371)*V(21)
  A(372) = RCT(372)*V(87)
  A(373) = RCT(373)*V(176)
  A(374) = RCT(374)*V(26)*V(197)
  A(375) = RCT(375)*V(25)*V(197)
  A(376) = 1.7e-11*V(27)*V(197)
  A(377) = RCT(377)*V(105)*V(189)
  A(378) = RCT(378)*V(105)*V(196)
  A(379) = RCT(379)*V(105)
  A(380) = RCT(380)*V(105)*V(189)
  A(381) = 2.3e-12*V(105)*V(195)
  A(382) = 2.63e-11*V(38)*V(197)
  A(383) = 3.78e-12*V(38)*V(195)
  A(384) = RCT(384)*V(116)*V(189)
  A(385) = RCT(385)*V(116)*V(196)
  A(386) = RCT(386)*V(116)
  A(387) = RCT(387)*V(116)
  A(388) = RCT(388)*V(116)*V(189)
  A(389) = RCT(389)*V(116)*V(189)
  A(390) = RCT(390)*V(116)*V(195)
  A(391) = RCT(391)*V(116)*V(195)
  A(392) = 4.65e-11*V(33)*V(197)
  A(393) = 1.25e-11*V(33)*V(195)
  A(394) = RCT(394)*V(135)*V(189)
  A(395) = RCT(395)*V(135)*V(196)
  A(396) = RCT(396)*V(135)
  A(397) = RCT(397)*V(135)*V(189)
  A(398) = 2.3e-12*V(135)*V(195)
  A(399) = RCT(399)*V(130)*V(189)
  A(400) = RCT(400)*V(130)*V(196)
  A(401) = RCT(401)*V(130)
  A(402) = RCT(402)*V(130)*V(189)
  A(403) = 2.3e-12*V(130)*V(195)
  A(404) = RCT(404)*V(136)*V(189)
  A(405) = RCT(405)*V(136)*V(196)
  A(406) = RCT(406)*V(136)
  A(407) = RCT(407)*V(136)*V(189)
  A(408) = 2.3e-12*V(136)*V(195)
  A(409) = 7.3e-11*V(72)*V(197)
  A(410) = RCT(410)*V(72)
  A(411) = 9.77e-11*V(60)*V(197)
  A(412) = RCT(412)*V(60)
  A(413) = RCT(413)*V(106)*V(197)
  A(414) = RCT(414)*V(106)*V(197)
  A(415) = RCT(415)*V(106)
  A(416) = RCT(416)*V(106)
  A(417) = 2.08e-12*V(39)*V(192)
  A(418) = 7.16e-11*V(69)*V(197)
  A(419) = RCT(419)*V(69)
  A(420) = 9.64e-11*V(64)*V(197)
  A(421) = RCT(421)*V(64)
  A(422) = RCT(422)*V(117)*V(197)
  A(423) = RCT(423)*V(117)*V(197)
  A(424) = RCT(424)*V(117)
  A(425) = RCT(425)*V(117)
  A(426) = 2.08e-12*V(34)*V(192)
  A(427) = 8.31e-11*V(99)*V(197)
  A(428) = RCT(428)*V(99)
  A(429) = 1.12e-10*V(89)*V(197)
  A(430) = RCT(430)*V(89)
  A(431) = RCT(431)*V(107)*V(197)
  A(432) = RCT(432)*V(107)*V(197)
  A(433) = RCT(433)*V(107)
  A(434) = RCT(434)*V(107)
  A(435) = 2.4e-13*V(45)*V(197)
  A(436) = RCT(436)*V(45)
  A(437) = 1.5e-12*V(37)*V(197)
  A(438) = 9e-13*V(70)*V(197)
  A(439) = 9e-14*V(70)*V(195)
  A(440) = RCT(440)*V(146)
  A(441) = RCT(441)*V(146)
  A(442) = RCT(442)*V(146)*V(189)
  A(443) = RCT(443)*V(146)*V(196)
  A(444) = RCT(444)*V(146)
  A(445) = RCT(445)*V(146)*V(189)
  A(446) = 2.3e-12*V(146)*V(195)
  A(447) = 1.53e-12*V(71)*V(197)
  A(448) = 3.13e-13*V(71)*V(195)
  A(449) = RCT(449)*V(150)*V(189)
  A(450) = RCT(450)*V(150)
  A(451) = RCT(451)*V(150)
  A(452) = RCT(452)*V(150)*V(196)
  A(453) = RCT(453)*V(150)*V(189)
  A(454) = RCT(454)*V(150)*V(189)
  A(455) = RCT(455)*V(150)*V(195)
  A(456) = RCT(456)*V(150)*V(195)
  A(457) = RCT(457)*V(140)*V(189)
  A(458) = RCT(458)*V(140)*V(196)
  A(459) = RCT(459)*V(140)
  A(460) = RCT(460)*V(140)*V(189)
  A(461) = 2.3e-12*V(140)*V(195)
  A(462) = RCT(462)*V(131)*V(189)
  A(463) = RCT(463)*V(131)*V(196)
  A(464) = RCT(464)*V(131)
  A(465) = RCT(465)*V(131)*V(189)
  A(466) = RCT(466)*V(131)*V(195)
  A(467) = RCT(467)*V(131)*V(195)
  A(468) = RCT(468)*V(131)*V(189)
  A(469) = RCT(469)*V(131)
  A(470) = 1.65e-11*V(111)*V(197)
  A(471) = RCT(471)*V(111)
  A(472) = 1.25e-11*V(75)*V(197)
  A(473) = RCT(473)*V(75)
  A(474) = RCT(474)*V(35)*V(197)
  A(475) = RCT(475)*V(35)*V(197)
  A(476) = 1.96e-11*V(132)*V(197)
  A(477) = RCT(477)*V(132)
  A(478) = RCT(478)*V(132)*V(195)
  A(479) = 7.3e-13*V(57)*V(197)
  A(480) = RCT(480)*V(57)
  A(481) = 1.89e-11*V(48)*V(197)
  A(482) = RCT(482)*V(48)
  A(483) = 9.2e-13*V(83)*V(197)
  A(484) = RCT(484)*V(83)
  A(485) = RCT(485)*V(83)
  A(486) = 3.57e-11*V(73)*V(197)
  A(487) = RCT(487)*V(73)
  A(488) = RCT(488)*V(73)
  A(489) = RCT(489)*V(88)*V(197)
  A(490) = RCT(490)*V(88)
  A(491) = 1.31e-12*V(44)*V(197)
  A(492) = 2.5e-11*V(68)*V(197)
  A(493) = RCT(493)*V(68)
  A(494) = 1.09e-12*V(49)*V(197)
  A(495) = 2.13e-11*V(59)*V(197)
  A(496) = RCT(496)*V(59)
  A(497) = RCT(497)*V(139)
  A(498) = RCT(498)*V(139)*V(192)
  A(499) = RCT(499)*V(139)*V(196)
  A(500) = RCT(500)*V(139)*V(189)
  A(501) = RCT(501)*V(139)*V(195)
  A(502) = RCT(502)*V(149)
  A(503) = RCT(503)*V(149)*V(189)
  A(504) = RCT(504)*V(149)*V(196)
  A(505) = RCT(505)*V(149)*V(189)
  A(506) = 2.3e-12*V(149)*V(195)
  A(507) = RCT(507)*V(124)*V(196)
  A(508) = RCT(508)*V(124)
  A(509) = RCT(509)*V(124)*V(189)
  A(510) = 2.3e-12*V(124)*V(195)
  A(511) = RCT(511)*V(74)
  A(512) = 1.27e-12*V(74)*V(197)
  A(513) = RCT(513)*V(46)
  A(514) = 4.42e-12*V(46)*V(197)
  A(515) = RCT(515)*V(125)*V(189)
  A(516) = RCT(516)*V(125)*V(196)
  A(517) = RCT(517)*V(125)
  A(518) = RCT(518)*V(125)*V(189)
  A(519) = 2.3e-12*V(125)*V(195)
  A(520) = RCT(520)*V(58)
  A(521) = 1.85e-12*V(58)*V(197)
  A(522) = RCT(522)*V(47)
  A(523) = 4.21e-11*V(47)*V(197)
  A(524) = RCT(524)*V(65)
  A(525) = 2.5e-11*V(65)*V(197)
  A(526) = 1.03e-11*V(61)*V(197)
  A(527) = 3.74e-11*V(85)*V(197)
  A(528) = RCT(528)*V(85)
  A(529) = RCT(529)*V(128)*V(189)
  A(530) = RCT(530)*V(128)*V(196)
  A(531) = RCT(531)*V(128)
  A(532) = RCT(532)*V(128)*V(189)
  A(533) = 2.3e-12*V(128)*V(195)
  A(534) = 1.87e-11*V(79)*V(197)
  A(535) = RCT(535)*V(79)
  A(536) = RCT(536)*V(138)*V(196)
  A(537) = RCT(537)*V(138)
  A(538) = RCT(538)*V(138)*V(189)
  A(539) = 2.3e-12*V(138)*V(195)
  A(540) = RCT(540)*V(67)
  A(541) = 3.2e-11*V(67)*V(197)
  A(542) = 1.38e-12*V(36)*V(197)
  A(543) = RCT(543)*V(36)
  A(544) = RCT(544)*V(118)*V(189)
  A(545) = RCT(545)*V(118)*V(196)
  A(546) = 2.3e-12*V(118)*V(195)
  A(547) = RCT(547)*V(118)*V(189)
  A(548) = RCT(548)*V(118)
  A(549) = 1.79e-12*V(43)*V(197)
  A(550) = 3.25e-11*V(62)*V(197)
  A(551) = RCT(551)*V(62)
  A(552) = 5.86e-12*V(86)*V(197)
  A(553) = RCT(553)*V(86)

! Aggregate function
  Vdot(1) = 0.781*A(161)-A(163)
  Vdot(2) = 0.007*A(199)+0.007*A(200)+0.007*A(201)+0.007*A(202)+0.007*A(204)-A(205)+0.005*A(353)+0.005*A(354)+0.005&
              &*A(355)
  Vdot(3) = A(54)+A(62)+A(67)+A(107)+A(116)+0.361*A(164)+0.5*A(175)+1.172*A(209)+0.858*A(210)+A(214)+2*A(217)+0.7*A(219)&
              &+2*A(220)+A(221)+A(222)+A(233)+A(239)+0.36*A(258)+0.63*A(259)+0.93*A(273)+0.93*A(274)+A(282)+0.5*A(284)&
              &+A(285)+A(286)+0.085*A(322)+0.4*A(369)
  Vdot(4) = 0.27*A(266)+0.21*A(268)
  Vdot(5) = 0.04*A(90)+0.04*A(91)+0.04*A(92)+0.04*A(93)+0.04*A(94)+0.04*A(95)+0.04*A(96)+0.04*A(97)+0.04*A(100)+0.04&
              &*A(101)+0.04*A(102)
  Vdot(6) = A(171)
  Vdot(7) = A(183)
  Vdot(8) = 0.35*A(271)
  Vdot(9) = 0.25*A(93)+0.15*A(156)+A(188)+A(189)+A(325)+A(326)+A(327)
  Vdot(10) = A(375)
  Vdot(11) = A(374)+A(376)
  Vdot(12) = 0.134*A(164)+A(379)+A(380)+A(381)+A(386)+A(388)+A(390)+A(396)+A(397)+A(398)+A(401)+A(402)+A(403)+A(406)&
               &+A(409)+A(410)+A(411)+A(412)+A(418)+A(419)+A(420)+A(421)
  Vdot(13) = A(51)+A(56)+A(57)
  Vdot(14) = A(61)+0.5*A(137)+A(174)+A(243)+A(287)+A(387)+A(389)+A(391)+A(407)+A(408)+A(427)+A(428)+A(429)+A(430)
  Vdot(15) = 0.53*A(374)
  Vdot(16) = 0.18*A(375)
  Vdot(17) = A(379)+A(380)+A(381)+A(386)+A(387)+A(388)+A(389)+A(390)+A(391)+A(396)+A(397)+A(398)+A(401)+A(402)+A(403)&
               &+A(406)+A(407)+A(408)+A(409)+A(410)+A(411)+A(412)+A(418)+A(419)+A(420)+A(421)+A(427)+A(428)+A(429)+A(430)
  Vdot(18) = A(525)+A(552)
  Vdot(19) = -A(6)-A(13)+A(29)
  Vdot(20) = -A(70)-A(71)
  Vdot(21) = A(328)-A(371)
  Vdot(22) = A(12)-A(38)-A(40)
  Vdot(23) = A(275)-A(276)
  Vdot(24) = A(109)-A(110)
  Vdot(25) = -A(375)
  Vdot(26) = -A(374)
  Vdot(27) = -A(376)
  Vdot(28) = A(20)-A(27)-A(35)
  Vdot(29) = 0.5*A(204)-A(208)+0.1*A(264)
  Vdot(30) = -A(42)+0.1*A(159)
  Vdot(31) = 0.578*A(202)+0.68*A(203)-A(206)+0.067*A(265)
  Vdot(32) = 0.272*A(202)+0.32*A(203)-A(207)+0.033*A(265)
  Vdot(33) = 0.18*A(375)-A(392)-A(393)
  Vdot(34) = A(392)+A(393)-A(426)
  Vdot(35) = A(440)-A(474)-A(475)
  Vdot(36) = A(541)-A(542)-A(543)
  Vdot(37) = A(414)+A(416)+A(423)+A(425)+A(432)+A(434)-A(437)
  Vdot(38) = 0.53*A(374)-A(382)-A(383)
  Vdot(39) = A(382)+A(383)-A(417)
  Vdot(40) = -A(74)-A(105)
  Vdot(41) = -A(68)-A(225)
  Vdot(42) = 0.41*A(155)-A(224)-A(362)
  Vdot(43) = A(544)-A(549)
  Vdot(44) = A(457)-A(491)
  Vdot(45) = -A(435)-A(436)+A(525)+A(552)
  Vdot(46) = A(499)-A(513)-A(514)
  Vdot(47) = A(504)-A(522)-A(523)
  Vdot(48) = A(443)-A(481)-A(482)
  Vdot(49) = A(462)-A(494)
  Vdot(50) = A(47)-A(48)-A(58)-A(59)
  Vdot(51) = -A(111)+A(113)+0.53*A(173)+0.5*A(185)+0.6*A(369)
  Vdot(52) = A(44)-A(49)-A(50)-A(55)
  Vdot(53) = A(152)-A(193)-A(347)
  Vdot(54) = A(153)-A(194)-A(348)
  Vdot(55) = A(129)-A(195)-A(349)
  Vdot(56) = -A(15)+A(57)+0.09*A(141)+A(337)+0.19*A(356)
  Vdot(57) = A(442)-A(479)-A(480)
  Vdot(58) = A(503)-A(520)-A(521)
  Vdot(59) = A(463)-A(495)-A(496)
  Vdot(60) = A(378)-A(411)-A(412)
  Vdot(61) = A(515)-A(526)+A(529)
  Vdot(62) = A(545)-A(550)-A(551)
  Vdot(63) = A(24)-A(25)-A(39)-A(41)
  Vdot(64) = A(385)-A(420)-A(421)
  Vdot(65) = A(507)-A(524)-A(525)
  Vdot(66) = 0.25*A(133)-A(148)+A(151)+0.09*A(259)
  Vdot(67) = A(536)-A(540)-A(541)
  Vdot(68) = A(458)-A(492)-A(493)
  Vdot(69) = A(384)-A(418)-A(419)
  Vdot(70) = A(417)-A(438)-A(439)
  Vdot(71) = A(426)-A(447)-A(448)
  Vdot(72) = A(377)-A(409)-A(410)
  Vdot(73) = A(452)-A(486)-A(487)-A(488)
  Vdot(74) = A(498)-A(511)-A(512)
  Vdot(75) = A(447)+A(448)-A(472)-A(473)
  Vdot(76) = A(65)-A(66)+0.19*A(161)-A(330)
  Vdot(77) = 0.75*A(93)+0.41*A(156)-A(197)-A(352)
  Vdot(78) = -A(1)-A(2)-A(3)-A(4)-A(5)+A(6)+A(30)+A(32)+A(34)
  Vdot(79) = A(521)+A(523)-A(534)-A(535)
  Vdot(80) = A(154)-A(198)-A(351)
  Vdot(81) = A(131)-A(192)-A(346)
  Vdot(82) = -A(16)+A(19)-A(31)+A(291)+A(297)+A(301)+A(305)+A(310)+A(314)+A(319)
  Vdot(83) = A(449)-A(483)-A(484)-A(485)
  Vdot(84) = A(121)-A(196)+0.85*A(280)+0.5*A(284)-A(350)
  Vdot(85) = A(516)-A(527)-A(528)+A(530)
  Vdot(86) = A(549)+A(550)-A(552)-A(553)
  Vdot(87) = 0.5*A(81)-A(322)-A(372)
  Vdot(88) = A(451)+A(453)+A(456)+A(483)+A(485)+A(486)+A(488)-A(489)-A(490)
  Vdot(89) = A(395)+A(400)+A(405)-A(429)-A(430)
  Vdot(90) = -A(60)-A(61)+0.85*A(123)+0.85*A(124)-A(368)
  Vdot(91) = 0.5*A(81)-A(323)-A(325)-A(327)-A(328)
  Vdot(92) = 0.272*A(125)-A(200)-A(203)-A(354)
  Vdot(93) = 0.1*A(126)-A(212)-A(358)
  Vdot(94) = 0.365*A(322)-A(324)-A(326)
  Vdot(95) = -A(211)+0.1*A(263)-A(356)
  Vdot(96) = -A(62)+0.134*A(164)+A(165)+0.5*A(175)+0.278*A(209)+0.142*A(210)+0.16*A(258)+0.057*A(266)+0.65*A(271)+0.07&
               &*A(273)+0.07*A(274)+A(292)
  Vdot(97) = A(127)-A(213)-A(214)+A(296)-A(359)
  Vdot(98) = 0.44*A(209)+0.473*A(210)-A(321)
  Vdot(99) = A(394)+A(399)+A(404)-A(427)-A(428)
  Vdot(100) = 0.17*A(234)+0.17*A(261)-A(281)-A(285)
  Vdot(101) = 0.654*A(211)+0.29*A(234)+0.29*A(261)-A(277)
  Vdot(102) = 0.628*A(125)-A(199)-A(202)-A(353)
  Vdot(103) = 0.037*A(125)-A(201)-A(204)-A(355)
  Vdot(104) = 0.5*A(267)+0.51*A(269)-A(270)-A(357)
  Vdot(105) = 0.47*A(374)-A(377)-A(378)-A(379)-A(380)-A(381)
  Vdot(106) = A(379)+A(380)+A(381)+A(387)+A(389)+A(391)+A(402)+A(403)+A(409)+A(410)+A(411)+A(412)-A(413)-A(414)-A(415)&
                &-A(416)
  Vdot(107) = A(396)+A(397)+A(398)+A(401)+A(406)-A(431)-A(432)-A(433)-A(434)
  Vdot(108) = 0.85*A(128)+0.41*A(157)+0.5*A(185)-A(215)-A(216)-A(360)
  Vdot(109) = 0.26*A(259)+0.21*A(266)+0.2*A(267)-A(286)-A(370)
  Vdot(110) = 0.044*A(180)-A(298)-A(299)-A(300)-A(301)
  Vdot(111) = A(438)+A(439)-A(470)-A(471)+A(517)+A(518)+A(519)+A(526)+A(527)
  Vdot(112) = 0.96*A(96)+0.96*A(101)+A(130)-A(190)-A(191)-A(345)
  Vdot(113) = A(112)-A(114)-A(160)-A(162)
  Vdot(114) = 0.055*A(226)-A(316)-A(317)-A(318)-A(319)-A(320)
  Vdot(115) = 0.32*A(75)+A(86)-A(118)+0.16*A(135)+0.75*A(145)+0.791*A(195)+0.32*A(230)+A(237)+A(253)-A(332)-A(333)&
                &+A(349)+0.32*A(361)
  Vdot(116) = 0.82*A(375)-A(384)-A(385)-A(386)-A(387)-A(388)-A(389)-A(390)-A(391)
  Vdot(117) = A(386)+A(388)+A(390)+A(407)+A(408)+A(418)+A(419)+A(420)+A(421)-A(422)-A(423)-A(424)-A(425)+A(427)+A(428)&
                &+A(429)+A(430)
  Vdot(118) = A(542)-A(544)-A(545)-A(546)-A(547)-A(548)
  Vdot(119) = 0.12*A(159)-A(311)-A(312)-A(313)-A(314)-A(315)
  Vdot(120) = 0.25*A(119)+0.25*A(120)+0.25*A(135)+0.25*A(136)+0.25*A(138)+0.24*A(139)+0.25*A(141)+0.25*A(142)+0.15&
                &*A(143)+0.25*A(145)+0.25*A(147)-A(149)+0.3*A(219)
  Vdot(121) = 0.022*A(182)+0.055*A(226)-A(302)-A(303)-A(304)-A(305)
  Vdot(122) = A(281)-A(282)-A(283)-A(284)
  Vdot(123) = 0.03*A(80)+0.09*A(140)+0.113*A(233)-A(278)-A(279)-A(280)
  Vdot(124) = A(489)-A(507)-A(508)-A(509)-A(510)
  Vdot(125) = A(502)+A(505)+A(506)-A(515)-A(516)-A(517)-A(518)-A(519)+A(522)
  Vdot(126) = A(272)-A(273)-A(274)-A(275)+A(276)
  Vdot(127) = 0.12*A(181)-A(306)-A(307)-A(308)-A(309)-A(310)
  Vdot(128) = A(520)-A(529)-A(530)-A(531)-A(532)-A(533)
  Vdot(129) = 0.96*A(92)+A(111)-A(113)-A(161)-A(369)
  Vdot(130) = 0.1*A(376)-A(399)-A(400)-A(401)-A(402)-A(403)
  Vdot(131) = A(437)-A(462)-A(463)-A(464)-A(465)-A(466)-A(467)-A(468)-A(469)
  Vdot(132) = A(441)+A(444)+A(445)+A(446)+A(473)+A(475)-A(476)-A(477)-A(478)+A(479)+A(480)+A(481)+A(482)+A(528)+A(531)&
                &+A(532)+A(533)
  Vdot(133) = 0.044*A(180)-A(293)-A(294)-A(295)-A(296)-A(297)
  Vdot(134) = 0.081*A(80)-A(258)-A(265)-A(365)
  Vdot(135) = 0.7*A(376)-A(394)-A(395)-A(396)-A(397)-A(398)
  Vdot(136) = 0.2*A(376)-A(404)-A(405)-A(406)-A(407)-A(408)
  Vdot(137) = -A(107)+0.15*A(155)+0.5*A(175)+A(178)+A(179)+A(218)+A(229)+A(242)+A(243)+A(244)+A(245)+A(246)+A(247)&
                &+A(248)+A(249)+A(250)+A(251)+A(252)+A(253)+A(255)+0.08*A(273)+0.08*A(274)+A(315)
  Vdot(138) = A(534)-A(536)-A(537)-A(538)-A(539)
  Vdot(139) = A(476)+A(478)-A(497)-A(498)-A(499)-A(500)-A(501)+A(511)+A(514)+A(537)+A(538)+A(539)
  Vdot(140) = A(435)-A(457)-A(458)-A(459)-A(460)-A(461)+A(474)
  Vdot(141) = 0.12*A(159)+0.11*A(180)+0.12*A(181)+0.326*A(182)-A(288)-A(289)-A(290)-A(291)-A(292)
  Vdot(142) = -A(171)-A(180)-A(183)
  Vdot(143) = -A(158)-A(159)-A(166)+0.102*A(180)+A(340)
  Vdot(144) = A(71)-A(72)+0.05*A(75)-A(119)+0.03*A(135)-A(153)-A(176)-A(178)+0.36*A(194)+0.05*A(230)+0.05*A(361)
  Vdot(145) = A(21)-A(28)-A(36)+A(52)+A(64)+0.6*A(89)+A(105)+A(117)+A(134)+A(169)+A(170)+0.5*A(185)+A(223)+A(225)+0.41&
                &*A(258)+0.3*A(259)+A(262)+A(279)+A(283)+A(383)+A(393)+A(439)+A(448)+A(478)
  Vdot(146) = A(422)+A(424)-A(440)-A(441)-A(442)-A(443)-A(444)-A(445)-A(446)+A(524)+A(535)
  Vdot(147) = 0.12*A(82)+0.21*A(258)+0.075*A(266)+0.06*A(267)+0.26*A(268)+0.16*A(269)-A(271)-A(364)
  Vdot(148) = 0.655*A(99)-A(100)-A(101)-A(102)
  Vdot(149) = A(472)-A(502)-A(503)-A(504)-A(505)-A(506)
  Vdot(150) = A(431)+A(433)-A(449)-A(450)-A(451)-A(452)-A(453)-A(454)-A(455)-A(456)+A(470)+A(540)
  Vdot(151) = A(206)+A(207)+A(208)-A(209)-A(210)
  Vdot(152) = A(45)-A(53)+0.25*A(97)+0.25*A(98)+0.25*A(102)+0.25*A(119)+0.25*A(120)+0.25*A(133)+0.25*A(135)+0.25*A(136)&
                &+0.5*A(137)+0.25*A(138)+0.24*A(139)+0.25*A(141)+0.25*A(142)+0.25*A(145)+0.25*A(146)+0.25*A(147)+0.24*A(259)&
                &+0.1*A(335)
  Vdot(153) = 0.4*A(89)-A(94)-A(95)-A(96)-A(97)
  Vdot(154) = A(70)+0.18*A(75)-A(86)-A(129)+0.09*A(135)-A(145)+0.209*A(195)+0.18*A(230)-A(237)-A(253)+0.18*A(361)
  Vdot(155) = A(94)+0.94*A(95)+0.7*A(97)+0.7*A(98)+0.919*A(100)+0.7*A(102)+0.134*A(164)-A(167)-A(169)+0.085*A(209)+0.088&
                &*A(210)+0.216*A(234)+0.06*A(259)+0.216*A(261)+0.48*A(278)+0.072*A(280)+0.085*A(322)+A(323)+A(324)-A(336)&
                &-A(337)-A(338)+A(371)+0.25*A(373)+A(379)+A(380)+A(381)+A(386)+A(388)+A(390)+A(396)+A(397)+A(398)+A(401)&
                &+A(402)+A(403)+A(406)+A(409)+A(410)+A(411)+A(412)+A(418)+A(419)+A(420)+A(421)
  Vdot(156) = -A(104)-A(131)-A(147)+A(166)+0.209*A(192)-A(241)-A(252)
  Vdot(157) = 0.9*A(264)-0.943*A(266)-A(267)+0.193*A(270)
  Vdot(158) = 0.15*A(83)+0.32*A(258)+0.019*A(266)+0.01*A(267)+0.44*A(268)+0.28*A(269)-A(272)-A(363)
  Vdot(159) = -A(73)-A(120)-A(154)+A(158)-A(177)-A(179)+0.209*A(198)
  Vdot(160) = 0.5*A(85)+0.96*A(90)+0.96*A(91)+0.94*A(94)+0.9*A(95)+0.66*A(97)+0.3*A(98)+0.894*A(100)+0.66*A(102)+0.075&
                &*A(128)+0.375*A(144)+0.791*A(192)+0.5*A(239)+0.42*A(258)+0.38*A(259)+0.3*A(266)+0.23*A(267)-A(287)-A(367)
  Vdot(161) = -A(85)-A(128)-A(144)+A(184)-A(239)-A(250)
  Vdot(162) = 0.9*A(265)-A(268)-A(269)+0.193*A(270)
  Vdot(163) = 0.6*A(89)-A(90)-A(91)-A(92)-A(93)+0.345*A(99)
  Vdot(164) = A(216)-A(217)-A(218)-A(219)-2*A(220)-A(221)-A(222)-A(223)
  Vdot(165) = 0.53*A(82)+0.578*A(126)+0.13*A(141)+0.36*A(142)-A(164)-A(165)+0.025*A(209)+0.088*A(210)+0.5*A(226)+0.234&
                &*A(234)+0.6*A(235)+0.24*A(259)+0.234*A(261)+0.9*A(263)+0.3*A(266)+0.23*A(267)+0.09*A(268)+0.06*A(269)+0.52&
                &*A(278)+0.078*A(280)+A(316)+A(317)+A(318)+A(319)-A(335)+0.26*A(356)+0.7*A(358)+A(364)+0.11*A(372)+0.5&
                &*A(373)+A(464)+A(465)+A(466)+A(494)+A(495)
  Vdot(166) = A(60)-A(78)+A(118)-A(123)-A(137)-A(231)-A(243)+A(471)+A(473)+A(543)
  Vdot(167) = -A(79)-A(124)+A(132)+A(134)-A(138)-A(232)-A(244)
  Vdot(168) = 0.009*A(80)+0.5*A(103)+0.25*A(146)-A(259)-A(264)-A(366)
  Vdot(169) = -A(82)-A(126)-A(142)+A(172)+0.054*A(180)+0.209*A(212)-A(235)-A(247)
  Vdot(170) = A(68)-A(69)+0.32*A(75)+A(115)-A(133)+0.16*A(135)-2*A(150)-2*A(151)-A(152)+0.44*A(156)+A(186)+0.36*A(193)&
                &+A(225)+0.32*A(230)-A(254)-A(255)+A(256)+A(331)+0.85*A(334)+0.32*A(361)+A(413)+A(415)+A(432)+A(434)+A(450)&
                &+A(454)+A(455)+A(477)+A(484)+A(487)+A(490)+A(497)+A(500)+A(501)+A(513)
  Vdot(171) = A(108)-A(109)+A(110)-A(115)+A(117)-A(156)-A(186)-A(188)+A(197)-A(256)+0.15*A(334)+A(342)
  Vdot(172) = -A(77)+A(106)-A(122)-A(136)-A(238)-A(249)+0.612*A(270)
  Vdot(173) = -A(141)+0.143*A(202)+0.493*A(204)+0.346*A(211)-A(234)-A(246)+A(260)-A(261)-A(262)-A(263)
  Vdot(174) = A(43)+A(45)+2*A(46)+A(50)-A(51)-A(52)+A(53)+A(55)-A(56)-A(57)+A(58)+0.05*A(63)+A(73)+0.57*A(77)+A(78)+0.82&
                &*A(80)+0.35*A(82)+0.122*A(83)+0.5*A(85)+0.07*A(87)+0.07*A(88)+0.75*A(97)+0.75*A(98)+0.75*A(102)+0.07*A(103)&
                &+A(104)+0.5*A(116)+0.75*A(119)+1.25*A(120)+0.15*A(123)+0.063*A(125)+0.102*A(126)+0.075*A(128)+0.75*A(133)&
                &+0.75*A(135)+0.95*A(136)+0.5*A(137)+0.75*A(138)+1.22*A(139)+0.75*A(140)+0.95*A(141)+0.89*A(142)+1.255&
                &*A(143)+0.625*A(144)+0.75*A(145)+0.785*A(146)+1.25*A(147)+0.59*A(157)+0.5*A(159)+0.19*A(161)+A(162)+0.732&
                &*A(164)+A(177)+0.456*A(180)+0.6*A(181)+0.12*A(182)+A(186)+2*A(187)+A(188)+A(189)+0.375*A(209)+0.375*A(210)&
                &+A(219)+0.445*A(226)+A(228)+A(229)+A(231)+0.747*A(233)+0.4*A(235)+0.143*A(236)+0.39*A(238)+0.5*A(239)+0.07&
                &*A(240)+A(241)+A(257)+0.62*A(258)+0.2*A(259)+0.15*A(266)+0.07*A(267)+0.7*A(268)+0.44*A(269)+0.35*A(271)&
                &+0.08*A(273)+0.08*A(274)+A(286)+A(288)+A(289)+A(290)+A(291)+A(321)+0.55*A(322)+0.9*A(335)+A(338)+A(341)&
                &+A(343)+A(344)+A(351)+0.71*A(353)+0.71*A(354)+0.71*A(355)+0.3*A(358)+0.5*A(359)+A(360)+0.5*A(363)+A(367)&
                &+A(368)+0.4*A(369)+A(370)+A(373)+A(459)+A(460)+A(461)+2*A(467)+2*A(468)+2*A(469)+A(493)+2*A(496)
  Vdot(175) = -A(112)+A(114)-A(116)-A(157)-A(187)-A(189)+A(215)-A(257)
  Vdot(176) = 0.7*A(87)+0.7*A(88)-A(89)-A(98)-A(99)+1.2*A(103)+0.6*A(146)+0.7*A(240)-A(373)
  Vdot(177) = -A(63)-A(64)+A(69)+A(73)+0.32*A(75)+0.86*A(77)+0.93*A(79)+A(104)+0.5*A(120)+0.15*A(124)+0.75*A(133)+0.16&
                &*A(135)+0.38*A(136)+0.5*A(138)+0.5*A(147)+A(148)+2*A(150)+A(151)+0.5*A(159)+A(177)+0.1*A(181)+0.64*A(193)&
                &+0.32*A(230)+A(232)+0.75*A(238)+A(241)+A(254)+A(255)+A(311)+A(312)+A(313)+A(314)-A(329)+A(347)+A(351)&
                &+A(352)+0.32*A(361)+A(450)+A(454)+A(455)+A(459)+A(460)+A(461)+A(471)+A(484)+A(487)+A(493)+A(508)+A(509)&
                &+A(510)+A(512)+A(528)+A(531)+A(532)+A(533)+A(537)+A(538)+A(539)+2*A(546)+2*A(547)+2*A(548)+2*A(551)+A(553)
  Vdot(178) = A(61)+0.35*A(82)+0.122*A(83)+0.5*A(85)+0.102*A(126)+0.075*A(128)+0.5*A(137)+0.29*A(141)+0.14*A(142)+0.255&
                &*A(143)+0.375*A(144)-A(168)-A(170)+A(174)+0.5*A(181)+0.88*A(182)+0.085*A(209)+0.088*A(210)+0.5*A(226)+0.234&
                &*A(234)+0.4*A(235)+0.143*A(236)+0.5*A(239)+A(243)+0.234*A(261)+0.9*A(263)+0.65*A(271)+0.07*A(273)+0.07&
                &*A(274)+0.52*A(278)+0.078*A(280)+A(287)+A(302)+A(303)+A(304)+A(305)+A(306)+A(307)+A(308)+A(310)+0.5*A(321)&
                &+0.55*A(322)+A(323)+A(324)-A(339)+0.58*A(356)+0.3*A(358)+A(363)+A(371)+0.5*A(372)+0.25*A(373)+A(387)+A(389)&
                &+A(391)+A(407)+A(408)+A(427)+A(428)+A(429)+A(430)
  Vdot(179) = -A(83)-A(84)-A(127)-A(143)+0.47*A(173)+0.5*A(185)+A(213)-A(236)-A(248)
  Vdot(180) = 0.058*A(80)+0.23*A(87)+0.23*A(88)+0.23*A(103)+0.3*A(139)+1.11*A(140)+0.115*A(146)+0.243*A(199)+0.513&
                &*A(200)+0.743*A(201)-A(226)+0.14*A(233)+0.23*A(240)-A(260)+0.275*A(353)+0.275*A(354)+0.275*A(355)+A(365)&
                &+A(366)
  Vdot(181) = -A(87)-A(88)-2*A(103)-A(130)-A(146)+A(183)+0.36*A(191)-A(240)-A(251)
  Vdot(182) = 0.728*A(83)+A(84)+0.16*A(120)+0.2*A(137)+0.18*A(141)+0.595*A(143)+A(160)+0.03*A(161)-A(174)-A(175)+0.65&
                &*A(179)+0.791*A(198)+0.085*A(209)+0.117*A(210)+A(214)+A(217)+0.7*A(219)+2*A(220)+A(221)+A(222)+0.216*A(234)&
                &+0.857*A(236)+0.42*A(259)+0.216*A(261)+0.27*A(266)+0.2*A(267)+0.09*A(268)+0.06*A(269)+0.85*A(273)+0.85&
                &*A(274)+0.5*A(277)+0.48*A(278)+0.072*A(280)+A(320)-A(344)+0.36*A(356)+A(359)+A(363)+0.39*A(372)+0.5*A(373)&
                &+A(491)+A(492)
  Vdot(183) = A(74)-A(75)-A(76)+A(105)-A(121)-A(135)+0.15*A(136)+0.209*A(196)-A(230)+0.3*A(238)-A(242)
  Vdot(184) = -A(80)-A(81)-A(125)-A(139)-2*A(140)+A(171)+0.75*A(199)+0.48*A(200)+0.25*A(201)-A(233)-A(245)
  Vdot(185) = A(76)+0.07*A(79)-A(106)+A(122)+0.25*A(144)+0.64*A(191)-A(361)
  Vdot(186) = 0.344*A(80)+0.035*A(87)+0.035*A(88)+0.035*A(103)+0.025*A(125)+0.18*A(139)+0.29*A(140)+0.018*A(146)-A(173)&
                &+0.3*A(180)-A(182)-A(184)-A(185)+0.294*A(233)+0.035*A(240)+A(293)+A(294)+A(295)+A(297)-A(343)+0.285*A(353)&
                &+0.285*A(354)+0.285*A(355)
  Vdot(187) = A(72)+0.14*A(75)+0.57*A(77)-A(108)-A(117)+0.75*A(119)+0.09*A(120)+0.033*A(126)+0.07*A(135)+0.54*A(136)&
                &+0.25*A(147)+A(149)+A(176)+A(178)+0.35*A(179)+0.64*A(194)+0.791*A(196)+0.13*A(230)+0.09*A(234)+0.57*A(238)&
                &+A(249)+A(250)+A(251)+A(252)+0.09*A(261)+A(282)+0.5*A(284)-A(331)+A(345)+A(346)+A(348)+A(350)+A(357)+0.13&
                &*A(361)
  Vdot(188) = 0.19*A(75)+0.187*A(126)-A(132)-A(134)+0.35*A(135)+0.25*A(138)+0.25*A(141)+0.25*A(142)+A(190)+0.791*A(212)&
                &+A(218)+0.19*A(230)+A(242)+A(244)+A(245)+A(246)+A(247)+A(248)+0.5*A(277)+A(285)-A(334)+0.19*A(361)
  Vdot(189) = -A(3)+A(4)-A(7)-2*A(9)-A(10)+A(11)-A(20)-A(23)+A(32)+A(33)+A(35)-A(43)-A(67)-A(69)-A(72)-A(73)-A(75)-A(76)&
                &-A(77)-A(78)-A(79)-A(80)-A(82)-A(83)-A(85)-A(86)-A(87)-A(91)-A(95)-A(100)-A(104)-A(115)-A(116)-A(210)&
                &-A(222)-A(223)-A(261)-A(262)-A(266)-A(268)-A(273)-A(278)-A(279)-A(282)-A(283)-A(289)-A(294)-A(299)-A(303)&
                &-A(307)-A(312)-A(317)-A(323)-A(324)-A(377)-A(380)-A(384)-A(388)-A(389)-A(394)-A(397)-A(399)-A(402)-A(404)&
                &-A(407)-A(442)-A(445)-A(449)-A(453)-A(454)-A(457)-A(460)-A(462)-A(465)-A(468)-A(500)-A(503)-A(505)-A(509)&
                &-A(515)-A(518)-A(529)-A(532)-A(538)-A(544)-A(547)
  Vdot(190) = A(42)-A(43)-A(44)-2*A(45)-2*A(46)-A(47)+A(48)+A(49)+A(59)+A(67)+0.19*A(75)-A(97)-A(102)+A(107)+0.5*A(116)&
                &-A(119)-A(120)-A(133)-0.91*A(135)-A(136)-A(137)-A(138)-A(139)-A(141)-A(142)-A(143)-A(144)-A(145)-A(146)&
                &-A(147)+0.44*A(155)+0.39*A(157)+0.28*A(159)+0.5*A(175)+A(176)+A(177)+0.051*A(180)-A(186)-A(187)-A(188)&
                &-A(189)+A(217)-A(219)+2*A(227)-A(229)+A(230)+A(231)+A(232)+A(233)+A(234)+A(235)+A(236)+A(237)+A(238)+A(239)&
                &+A(240)+A(241)+A(254)+A(256)+A(257)+A(309)+0.88*A(329)+0.3*A(330)+A(332)+2*A(333)+0.15*A(334)+A(342)+0.18&
                &*A(361)+A(362)+A(367)+A(423)+A(425)
  Vdot(191) = 0.95*A(63)+A(64)-A(65)+A(66)-A(67)+A(78)+0.93*A(79)+0.53*A(82)+0.5*A(116)+0.15*A(123)+0.15*A(124)+0.578&
                &*A(126)+0.3*A(137)+0.5*A(138)+0.36*A(142)-A(155)+A(162)+A(168)+A(170)-A(176)-A(177)-A(178)-A(179)+0.105&
                &*A(180)+0.28*A(181)+0.098*A(182)+A(187)-A(217)-A(218)+A(224)+0.445*A(226)-2*A(227)-A(228)-A(229)-A(230)&
                &-A(233)-A(234)-0.4*A(235)-A(236)-A(237)-A(238)-A(239)-A(240)-A(241)-A(242)-A(243)-A(244)-A(245)-A(246)&
                &-A(247)-A(248)-A(249)-A(250)-A(251)-A(252)-A(253)-A(254)-A(255)-A(256)+0.1*A(259)+A(309)+0.5*A(321)+0.085&
                &*A(322)+0.12*A(329)+0.7*A(330)+A(332)+0.85*A(334)+A(339)+A(341)+A(343)+A(344)+0.7*A(358)+A(360)+A(364)&
                &+A(368)+0.4*A(369)+0.11*A(372)+2*A(436)+A(490)+A(508)+A(509)+A(510)+A(535)+A(543)+A(553)
  Vdot(192) = A(3)-A(4)-A(5)+A(7)-A(8)+2*A(9)+2*A(10)-A(12)-A(21)+A(22)+A(23)-A(24)+A(25)+A(26)+A(27)-A(32)+A(34)+A(36)&
                &+2*A(37)+A(38)+0.95*A(39)+A(40)+A(41)+A(43)-A(47)+A(48)+A(59)-A(65)+A(66)+A(67)+A(69)+A(72)+A(73)+A(75)+2&
                &*A(77)+A(78)+0.93*A(79)+0.91*A(80)+0.88*A(82)+0.85*A(83)+1.5*A(85)+A(86)+1.3*A(87)+1.3*A(88)+A(90)+A(91)&
                &-A(92)+0.96*A(93)+2*A(94)+1.88*A(95)+0.7*A(97)+0.7*A(98)+0.3*A(103)+2*A(104)-A(109)+A(110)-A(111)-A(112)&
                &+A(113)+A(114)+A(115)+A(116)+0.925*A(128)+A(136)+0.375*A(144)+0.15*A(146)+A(147)+A(160)+A(190)+A(210)&
                &+A(222)+A(238)+0.5*A(239)+0.3*A(240)+A(241)+A(249)+A(250)+A(251)+A(252)+0.06*A(259)+A(261)+0.1*A(264)+0.1&
                &*A(265)+0.79*A(266)+0.88*A(268)+0.06*A(269)+1.85*A(273)+0.85*A(274)-A(275)+A(276)+A(278)+A(282)+A(286)&
                &+A(287)+A(289)-A(290)+A(294)-A(295)+A(299)-A(300)+A(303)-A(304)+A(307)-A(308)+A(312)-A(313)+A(317)-A(318)&
                &+A(323)+A(324)+0.7*A(330)+A(345)+A(346)+A(357)+A(361)+A(363)+A(364)+A(365)+A(366)+A(367)+0.6*A(369)+A(370)&
                &+A(373)+A(380)+A(381)+A(388)+A(389)+A(390)+A(391)+A(397)+A(398)+A(402)+A(403)+A(407)+A(408)+A(409)+A(410)&
                &-A(417)+A(418)+A(419)-A(426)+A(427)+A(428)+A(438)+A(439)+A(445)+A(446)+A(447)+A(448)+A(453)+A(454)+A(455)&
                &+A(456)+A(460)+A(461)+A(465)+A(466)+A(467)+A(468)+A(479)+A(480)+A(483)+A(484)+A(485)+A(491)+A(494)-A(498)&
                &+A(500)+A(501)+A(505)+A(506)+A(509)+A(510)+A(511)+A(512)+A(518)+A(519)+A(520)+A(521)+A(526)+A(532)+A(533)&
                &+A(538)+A(539)+A(546)+A(547)+A(549)
  Vdot(193) = A(1)-A(2)-A(7)-A(8)-A(14)-A(17)-A(29)-A(30)+0.25*A(93)-A(98)+0.15*A(155)+0.15*A(156)+0.15*A(157)-A(159)&
                &-A(162)-A(180)-A(181)-A(182)-A(226)-A(258)-A(259)-A(285)
  Vdot(194) = A(51)+A(52)-A(54)+A(56)+A(57)+0.05*A(63)+0.728*A(83)+A(84)+0.5*A(85)+A(90)+A(91)+0.5*A(116)+0.075*A(128)&
                &+0.29*A(140)+0.33*A(141)+0.595*A(143)+0.375*A(144)+0.39*A(157)+0.56*A(159)+A(160)+0.03*A(161)+0.505*A(164)&
                &+A(165)+2*A(167)+A(168)+2*A(169)+A(170)+0.5*A(175)+0.495*A(180)+0.56*A(181)+0.625*A(182)+0.5*A(185)+0.6&
                &*A(209)+0.678*A(210)+0.89*A(226)+0.09*A(234)+0.857*A(236)+0.5*A(239)+0.21*A(258)+0.09*A(261)+1.5*A(277)&
                &-A(288)-A(293)-A(298)-A(302)-A(306)-A(311)-A(316)+0.5*A(321)+0.55*A(322)+A(327)+0.88*A(329)+A(331)+A(333)&
                &+A(335)+2*A(336)+2*A(337)+A(338)+A(339)+A(340)+A(341)+A(343)+0.67*A(356)+0.5*A(359)+0.5*A(363)+A(367)&
                &+A(370)+1.89*A(372)+A(477)+A(512)
  Vdot(195) = A(5)+A(8)-A(10)-A(11)-A(12)-A(22)-A(26)+A(28)-A(33)-A(34)-2*A(37)+A(38)+0.05*A(39)+A(40)-A(52)+A(58)-A(64)&
                &-A(88)-A(89)-A(90)-A(94)-A(105)-A(117)-A(134)+0.81*A(161)+A(162)-A(166)-A(169)-A(170)-A(183)-A(184)-A(185)&
                &-A(225)+A(271)+0.15*A(273)+0.15*A(274)+A(290)+A(295)+A(300)+A(304)+A(308)+A(313)+A(318)+0.3*A(330)+0.4&
                &*A(369)-A(381)-A(383)-A(390)-A(391)-A(393)-A(398)-A(403)-A(408)-A(439)-A(446)-A(448)-A(455)-A(456)-A(461)&
                &-A(466)-A(467)-A(478)-A(501)-A(506)-A(510)-A(519)-A(533)-A(539)-A(546)
  Vdot(196) = A(14)+A(15)+A(16)-A(17)-A(18)-2*A(19)+A(22)-A(23)-A(24)-A(26)+0.95*A(39)+A(41)+A(43)-A(44)+2*A(46)+A(51)&
                &+A(52)+A(53)+A(54)+A(55)+2*A(56)+A(58)+A(62)+0.05*A(63)+A(69)+A(72)+A(73)+0.27*A(75)+0.82*A(80)+0.35*A(82)&
                &+0.85*A(83)+A(86)+0.8*A(87)+0.8*A(88)+A(90)+A(91)-A(93)-A(96)+0.5*A(97)+0.5*A(98)+0.934*A(100)-A(101)+1.2&
                &*A(102)+0.8*A(103)+A(119)+A(120)-A(121)-A(122)-A(123)-A(124)-0.937*A(125)-0.898*A(126)-A(127)-0.925*A(128)&
                &-A(129)-A(130)-A(131)+A(133)+0.64*A(135)+0.5*A(136)+0.3*A(137)+0.5*A(138)+1.1*A(139)+0.91*A(140)+A(141)&
                &+0.64*A(142)+1.7*A(143)+0.375*A(144)+A(145)+0.9*A(146)+0.5*A(147)+A(148)+A(149)+2*A(150)-A(152)-A(153)&
                &-A(154)-A(155)-A(156)-A(157)+0.28*A(159)+0.19*A(161)+0.773*A(164)+A(167)+A(169)+A(174)+A(176)+A(177)+0.237&
                &*A(180)+0.28*A(181)+0.141*A(182)+0.5*A(185)+A(186)+A(187)-0.175*A(209)+0.825*A(210)+0.7*A(219)-A(221)+0.445&
                &*A(226)+A(228)+0.27*A(230)+0.887*A(233)+A(234)+0.4*A(235)+0.85*A(236)+A(237)+0.5*A(239)+0.8*A(240)+A(254)&
                &+0.05*A(258)+A(261)-0.1*A(263)+0.61*A(266)-0.5*A(267)+0.69*A(268)-0.56*A(269)-A(274)+A(277)+A(278)-0.85&
                &*A(280)+A(282)-0.5*A(284)+A(285)+0.5*A(321)-A(325)-A(326)+A(329)+A(331)+1.73*A(335)+2*A(336)+A(339)+A(341)&
                &+A(343)+A(344)+A(345)+A(346)+A(347)+A(348)+A(349)+A(350)+A(351)+A(352)+0.985*A(353)+0.985*A(354)+0.985&
                &*A(355)+A(356)+A(357)+0.3*A(358)+A(359)+0.27*A(361)+A(363)+A(365)+A(366)+A(370)+0.89*A(372)+A(373)+0.53&
                &*A(374)+0.18*A(375)-A(378)+A(379)+A(380)+A(381)-A(385)+A(386)+A(387)+A(388)+A(389)+A(390)+A(391)-A(395)&
                &+A(396)+A(397)+A(398)-A(400)+A(401)+A(402)+A(403)-A(405)+A(406)+A(407)+A(408)+A(410)+A(412)+A(414)+A(415)+2&
                &*A(416)+A(419)+A(421)+A(424)+A(425)+A(428)+A(430)+A(433)+A(434)-A(443)+A(444)+A(445)+A(446)+A(451)-A(452)&
                &+A(453)+A(456)-A(458)+A(459)+A(460)+A(461)-A(463)+A(464)+A(465)+A(466)+A(467)+A(468)+A(469)+A(471)+A(473)&
                &+A(475)+A(477)+A(480)+A(482)+A(485)+A(488)+A(493)+A(496)-A(499)-A(504)-A(507)-A(516)+A(517)+A(518)+A(519)&
                &+A(528)-A(530)+A(531)+A(532)+A(533)-A(536)-A(545)+A(546)+A(547)+A(548)+A(551)+A(552)+A(553)
  Vdot(197) = 2*A(13)-A(14)-A(15)-A(16)+A(17)-A(18)-A(20)-A(21)-A(22)+A(23)-A(25)+A(26)-A(27)-A(28)+2*A(31)+A(35)+A(36)&
                &+0.05*A(39)-A(42)-A(49)-A(51)-A(53)-A(54)+A(55)-A(60)-A(62)-A(63)-A(68)-A(70)-A(71)-A(74)+A(84)+0.5*A(85)&
                &-A(99)-A(106)-A(107)-A(108)-A(118)+0.15*A(123)+0.15*A(124)+0.063*A(125)+0.68*A(126)+0.15*A(128)-A(132)&
                &-A(148)-A(149)+0.44*A(155)+0.44*A(156)+0.44*A(157)-A(158)+0.36*A(159)-A(160)-A(161)-0.773*A(164)-A(167)&
                &-A(168)-A(171)-A(172)-A(173)-A(174)+0.243*A(180)+0.36*A(181)+0.239*A(182)-0.36*A(191)-0.209*A(192)-0.36&
                &*A(193)-0.36*A(194)-0.209*A(195)-0.209*A(196)-A(197)-0.209*A(198)-0.875*A(199)-0.74*A(200)-0.625*A(201)&
                &-0.15*A(202)-0.5*A(204)-A(206)-A(207)-A(208)+1.125*A(209)+0.125*A(210)-0.346*A(211)-0.209*A(212)-A(213)&
                &-A(215)-A(216)+A(217)+0.7*A(219)+2*A(220)+2*A(221)+A(222)-A(224)+0.89*A(226)+0.05*A(258)+0.37*A(259)-A(260)&
                &+0.9*A(263)-A(264)-A(265)+0.5*A(267)+0.5*A(269)-0.388*A(270)-A(271)-A(272)+A(274)-A(277)+0.15*A(280)-A(281)&
                &+0.5*A(284)+A(285)-A(286)-A(287)-A(321)-0.365*A(322)+A(323)+A(324)+A(327)+0.07*A(335)+A(345)+A(346)+A(347)&
                &+A(348)+A(349)+A(350)+A(351)+A(352)+0.985*A(353)+0.985*A(354)+0.985*A(355)+A(356)+A(357)+A(358)+A(359)&
                &+A(360)+A(362)+A(368)+2*A(371)+2*A(372)-A(374)-A(375)-A(376)-A(382)-A(392)-A(409)+A(412)-A(413)-A(414)&
                &-A(418)+A(421)-A(422)-A(423)-A(427)+A(430)-A(431)-A(432)-A(435)-A(437)-A(438)-A(447)-A(470)-A(472)-A(474)&
                &-A(475)-A(476)-A(479)+A(482)-A(483)+A(487)+A(488)-A(489)-A(491)+A(493)-A(494)+A(496)-A(512)+A(513)-A(514)&
                &-A(521)+A(522)+A(524)-A(526)+A(528)-A(534)+A(540)-A(542)-A(549)+A(551)-A(552)
  Vdot(198) = 0.476*A(80)+0.035*A(87)+0.035*A(88)+0.035*A(103)+0.038*A(125)+0.28*A(139)+0.45*A(140)+0.018*A(146)-A(172)&
                &+0.2*A(180)-A(181)+0.453*A(233)+0.035*A(240)+A(298)+A(299)+A(300)+A(301)-A(340)-A(341)-A(342)+0.425*A(353)&
                &+0.425*A(354)+0.425*A(355)
      
END SUBROUTINE Fun

! End of Fun function
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~



END MODULE aromatics_kpp_Function

